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The concept of fragility provides a possibility to rank different supercooled liquids on the basis 
of the temperature dependence of dynamic and/or thermodynamic quantities. We recall here the 
definitions of kinetic and thermodynamic fragility proposed in the last years and discuss their inter- 
relations. At the same time we analyze some recently introduced models for the statistical properties 
of the potential energy landscape. Building on the Adam-Gibbs relation, which connects structural 
relaxation times to configurational entropy, we analyze the relation between statistical properties 
of the landscape and fragility. We call attention to the fact that the knowledge of number, energy 
depth and shape of the basins of the potential energy landscape may not be sufficient for predicting 
fragility. Finally, we discuss two different possibilities for generating strong behavior. 
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INTRODUCTION 

Soon after the introduction of the concept of "topo- 

Kaphic view of the Potential Energy Landscape (PEL)" 
, 2J , it became immediately clear that a key role in con- 
troUing the kinetic arrest characterizing the glass tran- 
sition was played by the number of distinct PEL 
local minima (inherent structures), 51 at, and by their 
energy distribution, ^^{E). Indeed, it was suggested 
that the qualitatively different behavior of different su- 
percooled liquids could be traced back to the differ- 
ence in the il,N function, or, more specifically to the 
steepness of the N dependence of this quantity. From 
general arguments, in a mono-component collection of 
a large number, TV, of units (atoms, molecules, ...), it 
can be shown that Ojv ~ exp(aA^). Similarly, it holds 
V,]\[{E) ~ exp (E(e)/fcB). Here E(e) assumes the mean- 
ing of "configurational entropy" and it is an extensive 
function of the energy per particle e = E/N. The quan- 
tity a {a — maxei^ie)} /NUb) is a measure of the total 
number of " inherent structure" (individual minima of the 
potential energy hyper-surface) . In comparing the behav- 
ior of different glass forming systems, particular empha- 
sis is placed in the relation existing between a and the 
" fragility" of the system under investigation. 

The "fragility" concept, in its modern form, has been 
introduced, developed and widespread by Angell 0]. It 
describes, in its kinetic version, how fast the structural 
relaxation time (tq,) increases with decreasing temper- 
ature on approaching the glass transition temperature, 
Tg, defined as the temperature where Tq, becomes equal 
to 100 s. "Strong" systems (low values of fragility) show a 
"weak" T-dependence of Tq,(T), which can be described 
by an Arrhenius law (Ta{T) = Too exp (A/fc^T)) while 
" Fragile" systems show -close to Tg- a much faster T de- 
pendence of the relaxation time, which is also markedly 
non-Arrhenius (this dependence could be, for example, 
described by a T- dependence of the activation energy 
A). The relaxation time is a quantity which is rather dif- 
ficult to access, in particular when the value of Tq is large, 



and, moreover, it seems also to be technique-dependent. 
For these reasons, in non-polymeric liquids, the fragility 
is usually defined through the T dependence of the shear 
viscosity, rj |^. This choice leads to a first ambiguity, 
especially in comparing different systems, as the fragility 
defined through Tq,(T) and that defined through 77 (T) are 
not coincident. This can be rationalized by recalling the 
Maxwell relation, 77=6*00 Tq (here Goo is the infinite fre- 
quency shear modulus of the liquid), and recalling that 
Goo at Tg spans over about two decades among different 
systems. Another possible definition of fragility comes 
from the temperature dependence of the mass diffusion 
coefficient. In this case, according to the Stokcs-Einstein 
relation {D = ksT/ (GTrrrj), being r the effective hydro- 
dynamic radius), it is the mobility fi {—D/T) that is (in- 
versely) proportional to the viscosity and, therefore, must 
be analyzed. Once more, it should be expected that the 
fragility defined via mobility and that defined via viscos- 
ity are not coincident. Indeed, i) the effective hydrody- 
namic radius may have a temperature dependence and 

ii) it is well known that in supercooled liquid at low tem- 
perature the " decoupling" phenomenon (the failure of the 
Stokes Einstein relation) occurs. In the recent years, the 
fragility has been quantified according to the T behavior 
of 77, but this has been done following different prescrip- 
tion {vide infra). 

Despite minor ambiguities introduced by its different 
definitions, the concept of fragility has a deep influence on 
the study of relaxation processes in supercooled liquids. 
Many studies have evidenced the existence of correlations 
between the values of the fragility and other properties 
of the supercooled liquids, such as: i) the "visibility" of 
the Boson Peak 0, Q ; ii) the T-dependence of the shear 
elastic modulus in liquids (shoving model) 0, 0, 0, 0| ; 

iii) the stretching of the decay of the correlation func- 
tions at the glass transition temperatur e [l^ ; iv) the 
nonlinearity of the relaxation functions |l4| and, very re- 
cently, v) the vibrational properties of the glass at T ^ 
[Tsf . Other works have tried to extract physical infor- 
mation on the nature of the glass transition from the 
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existence of these correlations 0|. Finally, we re- 
call a recent attempt to extend the dimensionality of the 
space spanned by the fragility index. Instead of using a 
single value to classify the T-dependence of the viscos- 
ity, Ferrer et al. 0| proposed to associate two indexes 
to every glassformer. The first index (fragile/non- fragile) 
measures how much the viscosity is Arrhenius- like at low 
temperature while the second one (strong/weak) does the 
same around the melting point. A deeper discussion on 
the correlation between fragility and other supercooled 
liquid properties can be found in Ref . [lij . 

The relation between the statistical properties of the 
landscape and the fragility is thought to be a central is- 
sue in the comprehension of the physics behind the glass 
transition. Debenedetti and Stillinger state in a very 
recent review: "Equally important is the translation of 
qualitative pictures ... into precise measures of strength 
and fragility based on the basin enumeration function" . 
A first connection between the fragility and the topo- 
graphic differences in the energy landscape is found in 
Ref. [21]. There the landscapes of strong liquids were 
supposed to have a "uniform" roughness, while a two- 
lengthscale arrangement of the minima -with the intro- 
duction of the meta-basins, a concept that has been re- 
cently revitalized by Doliwa and Heuer was expected 
to characterize the PEL of fragile liquids. In 1995, Angell 
|2,ll |. rationalizing the much larger specific heat jump at 
the glass transition shown by the fragile liquids with re- 
spect to the strong ones, concluded that "Fragile liquids 
would have high density of minima per unit energy..." 
and "Surfaces with few minima ... generate str ong liq- 
uids...". Similar conclusions are rep orted in Ref. |24| and 
by Debenedetti and Stillinger [20| who, more recently, 
wrote that "... strong landscape may consist of a single 
metabasin whereas fragile ones display a proliferation of 
well-separated metabasin" . 

Summing up, there seems to be consensus on the state- 
ments 



strong systems 
fragile system 



small a 
large a 



An attempt to determine a quantitative relation be- 
tween fragility and number of states on a theoretical 
basis, within the framework of the "gaussian landscape 
model" (see below), is due to Speedy 0, who derived 
a direct proportionality between kinetic fragility and a. 
This relation has been then criticized by Sastry [2g , who 
-again using the gaussian model to fit his molecular dy- 
namics simulation of the Kob-Anderson Lennard Jones 
Binary Mixture (BMLJ) at different densities- reached 
the conclusion of a proportionality between fragility and 
the square root of a. 

In this paper we first present a summary of the differ- 
ent definitions of "fragility" that are commonly used in 
the current literature, and then recall several models of 



configurational entropy (several "landscapes") proposed 
in the past that -with the help of the Adam-Gibbs equa- 
tion, or of the Vogel-Tamman-Fulcher relation, or both- 
lead to different expression for the fragility in terms of 
the parameters characterizing the "landscapes". In the 
subsequent sections, we review the Speedy and the Sas- 
try propositions on the a- dependence of the fragility 
for the examined landscapes. Finally we emphasize that 
landscapes with the same statistical properties (i.e. same 
total number of basins, same energy distribution of the 
basins depth) may be characterized by different fragili- 
ties, calling attention on the role of the different param- 
eters entering in the Adam- Gibbs expression. We con- 
clude discussing the obtained results in the context of the 
strong-to-fragile transition observed in some strong glass 
forming liquids. 



FRAGILITIES 

As discussed in the introduction, and following An- 
gell yi, we will define the kinetic fragility in terms of 
the temperature behavior of the viscosity and not of the 
structural relaxation time. Having clarified this point, 
however, we have to face -for the present purpose- differ- 
ent definitions of the " index of (kinetic) fragility" . The 
robustness of a concept like the fragility lies in the ob- 
servation that -when plotting log(?7(r)) vs. T/Tg- the 
curves for different liquids (beside very few exceptions) 
do not intersect each other, and converge to a common 
point dX T = Tg (by definition) and at T ^ cxi. Given 
this situation, it is possible to sort the systems, i. e. to 
unambiguously asses whether or not a system is more 
fragile than another. It is, therefore, natural to assign a 
numerical value to this concept: the index of fragility. 



Kinetic fragility: local definitions 

The first definition, let's call it "Angell's kinetic 
fragility" , , is 



dlog(r7(T)/77oo) 



d{Tg/T) 



T=T„ 



(1) 



Here 7700 is the limiting high temperature viscosity and Tg 
is defined from the condition rj{Tg) — 10^^ poise. As it is 
experimentally observed that all the liquids share a very 
similar value of 7700 — lO"** poise, this quantity is conven- 
tionally fixed to this value. Accordingly, an ideal strong 
glass (strictly Arrhenius behavior) would have = 17, 
whereas higher values are indication of higher fragility. 
While in principle there is no upper limit for m^, on a 
practical ground the most fragile systems seems to be 
tri-phenyl-phosphate, with « 160. 
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A very similar definition has been proposed by Speedy 



log(»)(T)/»)„ 



log('j(Ts)/»)o 



d{TjT) 



(2) 



At a first sight, it seems that a trivial normalization fac- 
tor would bring from to . However, this expression 
become more useful than Eq. 1 if we want to relax the 
assumption rjoa — 10^** poise. In conjunction with Eq.|21 
it is also useful to define the glass transition temperature 
Tg as the temperature where 'r]{Tg)/r]oo = 10^^; we will 
use this definition hereafter. As we will see below, if we 
aim to study, for example, the density dependence of the 
fragility of a given system, it will be easier to use Eq. 2 
where the density dependence of rjoo , although small, has 
been washed out. It is worth to point out, however, that 
for all the practical purposes, when dealing with the ex- 
perimental data the difference in using Eq. 1 or Eq. 2 is 
by all means irrelevant (apart from a trivial factor very 
close to 17). The fragility index ranges from one for 
strong glasses to wlO for the more fragile systems. 

The previous two definitions focus on the behavior of 
r]{T) at the glass transition temperature. More recently, 
another index of fragility — often referred to as F1/2 — 
has been introduced by Richert and Angell to "mea- 
sure" the fragility at intermediate temperature (see also 
the discussion in Ref. 13). Naming T* the temperature 
that satisfies log(77(T*)) = [log(r;(Tg)) + \og{r]iT^))]/2 
(i. e. the temperature where the viscosity is halfway -in 
logarithmic scale- between r]oo and 10"^^ poise), is de- 
fined as Fi/2 = 2(Tg/T*) — 1. It is worth to mention that 
Fi/2 and (or m^) do not provide "exactly" the same 
information: a plot of one quantity against the other 
does not indicate a perfect correlation, rather it shows a 
scatter of the points around an average trend H^. The 
existence of such a scattering has been recently rational- 
ized by Chandler and Garrahan within the framework of 
a coarse-grained model of glass formers [sO] • 

Finally, a generalized, temperature dependent fragility 
(either or m^) is sometime introduced, using equa- 
tions similar to Eq.s or |21 where Tg is substituted by 
a generic reference temperature T. We will call these 
quantities as 'm^{T) and mg{T), with the implicit defi- 
nition that when the argument is missing, the quantities 
are calculated at T = T^. 



Kinetic fragility: global definitions 

The previous indexes of fragility were associated to the 
behavior of r]{T) at a given temperature. Other defini- 
tions are based on the global behavior of the viscosity, 
and necessarily rely on the existence of a functional ex- 
pression for ri{T). 



A global definition of kinetic fragility arises from the 
experimental observation that the temperature depen- 
dence of the viscosity follows rather closely a Vogel- 
Tamman-Fulcher (VTF) law 'sH: 



viT) = Voo exp( 



T-Tr, 



(3) 



where ?7oo, D and To are system dependent parameters. 
As long as the VTF description of T?(r) is correct, one 
of the two parameters in the argument of the exponen- 
tial can be eliminated in favor of Tg as -from the defini- 
tion of glass transition temperature- the following rela- 
tion holds la: 



Tg = nil 



D 



17 In(lO)' 



(4) 



Plugging Eq.|21in Eq.[21 and using Eq.0I one gets that the 
parameter D is related to the previously defined fragili- 
ties: 



D = 



17 In(lO) 



(5) 



and, therefore, can be assumed to be a further fragility 
index. This index, which ranges from 00 for strong liquids 
(actually D « 100 for vitreous silica) to « 5 for the 
fragile ones, is in same sense "weaker" than the other 
three previously introduced, as its validity is based on 
the assumed T-dependence of the viscosity (EqlSJ. 

The assumption of the validity of the VTF law for 
the viscosity also leads to a relation between the local 
fragility defined at different temperatures. Indeed, re- 
calling the definition of F1/2 and Eq. [SJ one gets [23 



/2 



(6) 



Thermodynamic fragility 

An important step forward in relating the fragility with 
the PEL properties has been certainly achieved with the 
introduction of the "thermodynamic fragility" [32| . Simi- 
larly to the kinetic fragility which naturally emerges from 
the Angell plot (log(r/) vs. Tg/T for different systems), 
the vigor of the concept of thermodynamic fragility arises 
from the temperature dependence of the excess entropy 
S^^{T), defined as the difference between the entropy of 
the liquid and the entropy of the stable crystal. On plot- 
ting S^^{Tg)/S^^{T) vs. Tg/T, one obtains a plot very 
similar to the Angell plot, where the different systems 
stand in the same order • 

In similar fashion to the kinetic fragility it has 

been defined a "thermodynamic" fragility -F3/4: naming 
T* the temperature where S^^(Tg)/ S^^{T) = 3/4, i. e. 
the temperature where the inverse excess entropy equals 
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3/4 of its Tg value, F3/4 is defined as F3/4 = 2{Tg/T*) - 
1. In this case, the value 3/4, and not 1/2, has been 
chosen because of the difficulties associated to determine 
the excess entropy at high T/Tg in strong liquids. In a 
recent paper Martinez and Angell have shown that 
it exists a remarkable correlation between F1/2 and -F3/4: 
with few exceptions it turns out that ~ F3/4 within 
10%. This observation rationalizes the well known fact 
that the amplitude of the specific heat jump at Tg is 
linked to the fragility, but also points out that is not the 
specific heat jump alone, but rather this jump divided by 
the excess entropy at Tg , that is actually related to . 

In analogy with (or with ) it would be natural to 
define a further index of the thermodynamic fragility as 
the derivative at Tg of the inverse reduced excess entropy 
with respect to the inverse reduced temperature. To our 
knowledge, this index has not been yet introduced, but 
-as we will see below- this quantity naturally appears 
when the Adam-Gibbs relation is used to work out a link 
between kinetic and thermodynamic fragility. It is use- 
ful, therefore, to introduce this thermodynamic fragility 
(m^) index as: 

being 5'^ (T) the temperature derivative of S^^ (T) . 



and Vogel-Tamman-Fulcher relations (Eq. ^ are valid, 
then necessarily To and Tx are equal one to each other. 
This equality has been recently disputed We do 

not further discuss this problem, with the aim to study 
the mathematical consequences of the different landscape 
models introduced in the literature, we will assume (when 
necessary) that f is a slowly varying smooth function of 
T (thus, that To = Tk)- It must also be noted that 
the thermodynamic fragility is defined through the ex- 
perimentally accessible excess entropy, while the Adam- 
Gibbs relation calls into play the configurational entropy. 
In the following we will not make difference between the 
two entropies, relying upon the observation that config- 
urational and excess entropy seems to be actually pro- 
portional to each other even if other studies indi- 
cate the failure of such a proportionality . Assuming 
that the Adam-Gibbs relation correctly describes the T- 
dependence of the viscosity in a supercooled liquids, by 
plugging Eq. iniinto the definition of m^, Eq. |21 we get 
(using viTg)/Voo = IQi^): 

and, recalling Eq.Q we have the desired relation between 
kinetic and thermodynamic fragility: 



Relation between 
kinetic and thermodynamic fragility 

The Adam-Gibbs equation [s^l establishes a relation 
between the structural relaxation time and the configu- 
rational entropy S(T): 

r(T)=ro,exp(^;^) (8) 

or, relying on the Maxwell relation, between the viscosity 
and the configurational entropy: 

77(T)=7yo,CXp(^;^) (9) 

where t^o ('yoo) is the usual infinite temperature limit for 
the relaxation time (viscosity) and £ a system dependent 
parameter with the physical dimension of an energy that 
is somehow related to the energy barrier for activated 
processes. This equation is the key relation that allows 
us to create a link between kinetic and thermodynamic 
fragility and, ultimately, via the configurational entropy 
a link between kinetic fragility and the statistical prop- 
erties of the PEL. Let us first observe that, as the energy 
barrier is expected to have a weak and smooth temper- 
ature behavior and not to diverge at any temperature, 
according to Eq. El the viscosity diverges at the temper- 
ature (Kauzmann temperature Tk) where the configura- 
tional entropy vanishes. If both the Adam-Gibbs (Eq. O 



TOg = 1 + TOy . (11) 

EauationllOlalso constitutes the basis to obtain a link be- 
tween the kinetic fragility and the number of states a. 
Indeed, recalling the relation a = maXe{^{e{T))}/NkB, 
if we know -or have a model for- the configurational en- 
tropy of a given system, we could determine a and , 
and thus try to relate one to the other. 



MODELS OF LANDSCAPE 

In this section we will briefly recall the main mod- 
els that have been introduced in the recent literature to 
represent the configurational entropy of supercooled liq- 
uid systems. In the first three subsections we elucidate 
models of configurational entropy and derive the relations 
between the different quantities of interest (T and e de- 
pendence of S, fragility, etc.) with the specific hypothesis 
that the vibrational entropy associated to a specific mini- 
mum of the PEL is independent from its energy elevation. 
In the following subsection, we relax this hypothesis, as- 
suming a linear dependence of the vibrational free energy 
from e, and showing how the equations relating the rel- 
evant physical quantities to the configurational entropy 
parameters are modified. 
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Gaussian model 

The gaussian model is at the basis of the interpreta- 
tion of the configurational entropy in simulated super- 
cooled hquids. After the first studies [s^ Is^ 0|, the 
gaussian model has been chosen to describe quantita- 
tively the ener gy d ependence of E(e) in different systems 
HE HE ^ lialia l . According to this model, an explicit 
functional form (gaussian) for ^^{E) -the energy distri- 
bution of the minima of the PEL- is assumed, 



riNiE) = exp(aA^) exp[- 



{E - Eo) 



(12) 



From this equation, the configurational entropy of the 
gaussian model becomes (e = E/N): 



S(e) = UbN 



(13) 



being e = e/VN. In this expressions a counts the total 
number of states (it is the maximum of Y,{e)/N in ks 
units), Co is an irrelevant parameter (it fixes the zero of 
the energy scale) and e is the width of the distribution. In 
order to express the configurational entropy as a function 
of the temperature, we must first determine the energy of 
the minima of the PEL populated at a given temperature. 
Using lillil 



1 _ dS(e)/iV 
T ~ de 



we get 



e(T) 



2kBT 



(14) 



(15) 



and, finally, inserting Ea. 1151 into Eg. 1131 we have the ex- 
plicit expression of the configurational entropy as a func- 
tion of the temperature: 



I](T) - keN 



{2kBTY 



(16) 



From Eg. 1131 the Kauzmann energy ex, i-C. the energy 
where S(e) = 0, is promptly derived: 

(17) 



eK = e-o - ^\IOi 



and, plugging the Kauzmann energy (Eq. [T7|) in Eq. 
we find the Kauzmann temperature: 



^bTk = 



2^' 



(18) 



It is useful to eliminate e from the expression of the 
configurational entropy (in its explicit T-dependent ex- 
pression) in favor of T^, using Ea. I18L so to obtain: 



I](r) = keNa 



rp2- 

2^2 



(19) 



Once we have a model for the configurational entropy, 
we can -applying Eq. llOt find an expression for the 
fragility in terms of the parameters of the model itself. 
As parameters, we have the freedom to choose among 
(a, e, Tff, ex)- One compact possibility, which has the 
advantage to explicitly depend only on Tk^ is: 



(20) 



In this expression, Tg appears explicitly and cannot be 
eliminated because in the gaussian model (a 'pure thermo- 
dynamic model) the dynamics is not defined and there- 
fore Tg must be regarded as a parameter external to the 
theory. Other possible selection of parameters, and thus 
other expressions for the fragility, are of course possible. 
Eq. 1201 (as well as similar expressions for other landscape 
models, see below) makes clear the well know fact that 
the fragility is somehow related to the " distance" between 
Tg and Tk'- the higher is the ratio Tg/Tx the strongest 
is the liquid. 

As finally remark, we observe how -having imposed the 
validity of both the Adam-Gibbs relation and the gaus- 
sian model for the configurational entropy- the temper- 
ature dependence of the viscosity turns out to be con- 
trolled by the law: 



V{T) = rjoo exp( 



DT 



K 



T 



T - Tk T + Tk 



(21) 



with D = £/{aNkBTK), which is different by a VTF re- 
lation. In other words, the VTF law, the Adam-Gibbs re- 
lation and the gaussian model cannot be simultaneously 
invoked (especially when the shape of the PEL basins is 
independent on the depth). Eguationl21lcan be regarded 
as a VTF law with a temperature dependent coefficient 
D'{T) = DT/{T + Tk). In the high T limit {T »Tk), 
D' ^ D while in the low T regime (T approaching T^) 
D' D/2. 



Hyperbolic model 

For thirty years it has been realized ^E| that the 
temperature dependence of the (constant volume) ex- 
cess specific heat can be described by a hyperbolic law 
(C ~ const + const' /T), and this law is commonly used 
to represent the experimental data [23|- The "landscape 
model" that gives rise to such a temperature dependence 
for the excess specific heat is the so called hyperbolic 
model, recently introduced and discussed in detail by 
Debenedetti, Stilhnger and Lewis [13 . In Ref. [i^l, the 
model is derived from the assumption of a hyperbolic 
temperature dependence of the "configurational" heat ca- 
pacity, and (assuming the validity of the Adam-Gibbs re- 
lation) , it implies as a mathematical consequence the va- 
lidity of the VTF relation. For simplicity, here we prefer 
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to start assuming the mathematical validity of both the 
Adam-Gibbs and the VTF, the hyperbolic temperature 
dependence of the excess specific heat results as conse- 
quence. Obviously, as discussed in 48] the two routes are 
equivalent. It is worth to point out that the "gaussian 
landscape" is named after the e-dependence of the num- 
ber of states, while the "hyperbolic landscape 
after the T behaviour of the specific heat |2^^ 
different quantity. It is our aim to write down the main 
expressions for this model using the same notation of the 
previous section, and to extract the equations for the 
fragilities. By comparing Eq. |31 and O it turns out an 
explicit temperature dependence for S(T): 



" is named 
a rather 



i](r) = 



E 



DTk 



1-^ 
T 



(22) 



It is implicit in this expression the coincidence of To and 
Tk- This equation can be cast in form very similar to 
Eg. 1161 bv defining the quantities a and e: 



E 



DNkeTK ' 



2£ 

Un 



= 2kBTKa 



(23) 



(24) 



As we will see soon, a and e play here the same role 
as they have in the gaussian model, therefore the first 
equation is a link between the "number of states" and 
the constants entering in the AG (£) and VTF {D and 
Tk) relations. The second equation can be compared to 
Eg. 1181 where ^/a appears instead of a. Rewriting Eg. 1221 
with the elimination of E and D in favor of a and e, we 
have: 



E(T) = keN 



2kBT 



(25) 



an expression that can be directly compared with Eg. 1151 
or, expressing the pre-factor in Eq. [S^in terms of a via 
Eq.Eni 



E(T) = ksNa 



1 - 



Tk 
T 



(26) 



that can be compared with Eq. [T^ 

At variance with the gaussian model, where we started 
with a model for S(e) and derived S(r), we now have a 
model for S(T). To obtain an expression for I](e) we first 
derive the temperature dependence of the energy of the 
minima visited by Eq. 1141 



e{T) 



T 



dT 



dT 



(27) 



cr + aksTK In (T/Tr) 



where en and Tr are integration constants whose values, 
as we will see, are not relevant for the interesting physical 



quantities. Inverting Eq. |27| and plugging the resulting 
r(e) into Eq.HHwe get: 



E(e) = Nksa 



1 



Tf 



K 



■ exp 



2(e-efl) 



(28) 



Obviously, we can eliminate Tr from this equation, by 
properly redefining cr. A useful possibility is to choose 
Tji = Tk, then, from Eq. |23 cr — ck and: 



I](e) = Nksa 



1 — exp 



2{e-eK) 



(29) 



At variance with the configurational entropy of the gaus- 
sian model, the present S(e) does not show any maxima, 
rather it increase continuously, asymptotically approach- 
ing the value Nksa. 

From Eq.[2Bl using Eq.^1 we can easily determine the 
fragihty of this model: 



T„ - Tk 



(30) 



It is worth to point out that this expression is the expan- 
sion of the fragility of the gaussian model to first order 
in Tg - Tk- 



Logarithmic (or binomial) model 

The previous two models for the configurational en- 
tropy share the property that dT,/de is non diverging 
at e = ck, so the Kauzmann temperature exists and 
it is non-vanishing. In order to introduce a more flexi- 
ble model, embedding the possibility of having a vanish- 
ing Kauzmann temperature, Debenedetti, Stillinger and 
Shell ^3 recently proposed a modification of the gaus- 
sian model that, with a slight change in notation with 
respect to the original definition, reads: 



E(e) = Nksa 
(1 



7 



{(1-7) 
ln(H- 



1 



) + (l 



)ln(l 



21n(2) 



with u < u < ^/a) given by: 



(31) 
(32) 



This is a linear combination -weighted by the param- 
eter 7- of the parabolic configurational entropy typical 
of the gaussian model and a term that depends on the 
logarithm of the energy. Here we want to describe in 
detail the properties of this model for the specific case 
7 = 1, i. .e. of a model that is totally "logarithmic". The 
logarithmic model is essentially a binomial distribution, 
i.e. model for the thermodynamics of a gas of binary 



7 



excitations |50j. It has been used to model the thermo- 
dynamics of supercooled liquids and the T-dependence 
of the inherent structure energy [s^l • Obviously, the log- 
arithmic term in Eq. 1311 become dominant in the low- 
T/low-(e— ck) region, therefore the model discussed in 
this section can be thought as an approximation of the 
Debenedetti, Stillinger and Shell model valid in the low 
T limit. It is, however, interesting to study such a model 
in the whole energy range. Indeed, as we will see below, 
a visual inspection of the function I](e) indicates that 
this model and the gaussian model represent very similar 
"landscapes", i. e. very similar distribution of the min- 
ima energy. Thus, we define the "logarithmic" landscape 
as: 



S(e)=iVfci3ax (33) 
(l + ^)ln(l + ^) + (l-^)ln(l-^)- 



2 ln(2) 



This expression for the configurational entropy has the 
properties to vanish at u = i\/a, therefore the Kauz- 
mann energy results to be at u = —\/a or explicitly 
GK = Bo — G\foi- At this energy, the derivative of I](e) 
shows a logarithmic divergence, thus implying that the 
Kauzmann temperature must vanish. Similarly to the 
gaussian model, the parameter Co is the energy of the 
"top of the landscape" and a represents the maximum 
of T,{e)/NkB- Using Eg. IH^ and the expression for ex, 
Eq. 1221 can be explicitly written in terms of the reduced 
energy measured with respect to the Kauzmann energy 
V = {e — gk)/^ as: 



E(e) = 7Vfcsa 



^ln(^) 



-(2- 



ln(2- 



21n(2) 



(34) 



We can now follow the same route used in the discussion 
of the gaussian model. Via Eq. ^1 with straightforward 
algebra, we obtain the temperature dependence of the 
energy of the minima: 



e(T) 



'atanh 



21n(2)e 



(35) 



and inserting this expression in Eq. 1341 the temperature 
dependence of the configurational entropy is promptly 
derived: 



2 cosh —p^' 



T 



e , / ln(2)e 

—= tanh J: ' 



} (36) 



As a consequence of the infinite value of 'S,{e)/de at 
ck, this function does vanish only at T = 0, i. e. for this 
model Tk = 0. It is convenient, for sake of compactness, 
to define a typical temperature which -in analogy with 
Tk in the gaussian and hyperbolic models- could be used 



to scale the temperatures in the logarithmic model. We 
arbitrarily introduce the quantity: 



T^ = iln(2)-^ 
3 ksva 



(37) 



whose value is very close to the "apparent" Kauzmann 
temperature that would have been identified by extrap- 
olating Eq. 1361 towards zero using only information on 
S(r) at "high" temperature, similarly to what is done 
experimentally. In other words, the logarithmic model 
predicts a temperature dependence of the configura- 
tional entropy that -around the inflection region- can 
be approximated by a straight line that goes to zero at 
kBTK^/a/e « 0.23 (« ln(2)/3). Having introduced the 
"apparent" Kauzmann temperature for the logarithmic 
model, we can write Eg .s 1351 and 1361 as : 



s(r): 



e(T). 



Nksa 
1^ 



afesT^tanhf ^ 



In 



ln(2) 



2 cosh 



(38) 



3r, 



K 



T 



tanh 



3T; 



K 



T 




0.5 1.0 1.5 

(e-ej^)/£ 



FIG. 1: Sketch of the energy dependence of the configura- 
tional energy for three models: gaussian (full line) , hyperbolic 
(dashed line) and logarithmic (dot-dashed line). The reduced 
entropy T.{e) /Nks is plotted as a function of (e — eA')/e for 
the specific case of q = 0.8 

Once the explicit T dependence of S(T) is known, both 
the fragility m ^ , definite in Eq.^J and the T dependence 
of the viscosity (from the Adam-Gibbs equation) can be 
worked out. The two expressions reads: 



/3Tj 



K 



\ T 



3T, 



K 



T 



cosh 



I cosh^ 
3T, 



3r, 



K 



K 



T 



T 



sinh 



In 
3T 



2 cosh 



m 



K 



T 



r 



(39) 
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gaussian model hyperbolic model 



logarithmic model 



E(e)/iVA:B{ 

ex 
ksTK 
e 

e{T) - Co 
e(T) - ex 

E{T)/NkB 



ln(77(r)/,7. 



Bo — 

2^ 

-2akBTK (^) 



1 - 



T-Tk T+Tk 



a[l — exp (—2m)] 



Q [1 — exp (—2d)] 



2akBTK 



aksTK In ( 



akuTxln { ^ 



a 1 



['-(*)] 



Tg-TK 



DTk 
T-Tk 



21n(2) 



_(1 )ln(l- 



[l-H^ [^M^) + (2-^)ln(2-^)]] 



Co — fL\fot 



ln(2) E 

3 "7^ 



T^)VakBT^ 



o / ST'* 

-j^akBTii tanh I ^ 



^afcsrl.[l-tanh(Hi:)] 



2cosh(^) -^tanhl^ 



■)} 



1+w 



'I cosh^ (ui) In [2 cosh(ui)] — ly cosh(u;) sinh(u;) | 



{MHli^^jln [2 



3 '-D^\ In 2cosh(^) - ^tanh 



TABLE I: Summary of the main relations relating the relevant quantities (left column) for the three configurational entropy 
models introduced before. The relation that define the model is reported in box. The variable u is defined as the reduced energy 
measured with respect to Co: u = (e — eo)/e, while v is that measured starting from the Kauzmann energy: v = (e — ex)/e. The 
variable w is a shortcut for 3T^/Tg. In the case of the logarithmic model, is used (see Eg. I37II in place of the Kauzmann 
temperature. 



rln(2) /3r^ 



3T, 



K 



To 



tanh ( ^ 



){'■ 



2cosh( 

T 



(40) 



Summary of models 

In Table I, we summarize the expressions derived in 
the framework of the three model examined before for 
different quantities. These quantities are: 

i) and ii) the configurational entropy as a function of e, 
in this case we explicitly report S(e) as a function 
of the variables u = (e — eo)/e and w = (e — e_R-)/eto 
emphasize that the zero of the energy is irrelevant 
and that e only acts as an energy scale. 



iii) The explicit expression of the Kauzmann energy in 

term of 60, a and e. 

iv) and v) The relations used to eliminate e in favor of 

Tk (or in the case of the logarithmic model). 

vi) and vii) The temperature dependence of the inherent 
structures energy, reported in terms of Tk ■ 

viii) The temperature dependence of the configurational 
entropy, now reported in term of Tk- 

ix) The expression for the fragility reported in terms of 

the thermal parameters. 

x) Finally, we report the temperature dependence of 

the viscosity resulting from the application of the 
model. In the last expression the parameter D is: 
D = £/{aNkBTK). 
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FIG. 2: Plot of the temperature dependence of the config- 
urational entropy for three models: gaussian (full line), hy- 
perbolic (dashed line) and logarithmic (dot-dashed line). The 
reduced entropy T.{T) /Nks further normalized to a is plotted 
as a function of T/Tk- In the case of the logarithmic model 
defined in Eg. 1371 is used in substitution of Tk- 




FIG. 3: Similarly to Fig.|21 the temperature dependence of the 
configurational entropy for three models (gaussian (full line), 
hyperbolic (dashed line) and logarithmic (dot-dashed line)) is 
plotted as a function of Tk/T. In the case of the logarithmic 
model Tj^- defined in Eg. 1371 is used in substitution of Tk. 



In Fig. n we sketched the e dependence of the con- 
figurational energy for the examined models: gaussian 
(full line) , hyperbolic (dashed line) and logarithmic (dot- 
dashed line). As an example, the three configurational 
entropies are reported for the specific case of a = 0.8 (as 
the scaling of S(e) with a for the hyperbolic model is dif- 
ferent from that for the gaussian and logarithmic models, 
we cannot use a reduced variable). 

Similarly, in Fig.s[21and|3|we report the corresponding 
configurational entropy as a function of T /Tk and Tk/T 
respectively. In the case of the logarithmic model 
defined in Eq. |^is used to scale the temperatures. 



In Fig. 0] we report the temperature dependence of 
the energy elevation (normalized to the factor akBTx) 
with respect to ck of the minima of the PEL visited 
at equilibrium for the three examined models: gaussian 
(full line) , hyperbolic (dashed line) and logarithmic (dot- 
dashed line). The hyperbolic model shows a non physical 
continued rise of e(T) at increasing T. 



2.5 




FIG. 4: Plot of the temperature dependence of the energy 
elevation for three models: gaussian (full line), hyperbolic 
(dashed line) and logarithmic (dot-dashed line). The energy 
elevation e — ck, normalized to ksTK, and further normalized 
to a, is plotted as a function of T/Tk- In the case of the 
logarithmic model T^ defined in Eq.|^is used in substitution 
of Tk. 

Finally, in Fig.|31we report in Arrhenius scale the tem- 
perature dependence of the viscosity for the three exam- 
ined models: gaussian (full line), hyperbolic (dashed line) 
and logarithmic (dot-dashed line). 



Gaussian models with 
non-constant vibrational entropy 

All the discussion in the previous sections was based on 
the assumption that the vibrational entropy associated to 
a given basin is independent from the energy elevation of 
the minimum of the basin itself. These assumption lead 
to the simplified microcanonical definition of tempera- 
ture reported in Eq. [141 Following recent experimental 
[stI and numerical (42j evidences indicating a vibra- 
tional entropy that actually depends on the energy of 
the minima, in the present section we relax the previ- 
ous assumption, and, for the specific case of the gaussian 
model, we develop the calculation in the case of an ex- 
plicit dependence of the vibrational entropy, S^, on e. 
In particular, taking advantage of the outcome of recent 
molecular dynamics calculations, we develop iS'i,(e) in se- 
ries of e — e/f and retain only the first order term, an 
approximation certainly valid for low enough tempera- 
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FIG. 5: Plot of the temperature dependence of the viscosity 
for three models: gaussian (full line), hyperbolic (dashed line) 
and logarithmic(dot-dashed line). The logarithm of the vis- 
cosity normalized by rjiToo), normalized to D is plotted as a 
function of Tk/T. In the case of the logarithmic model 
defined in Eg. 1371 is used in substitution of Tk- 



ture: 



Sv{g) — H ; — 

de 



(e - ex) 



(41) 



The quantity dSy/de is a further system-dependent pa- 
rameter. For sake of simplicity let us define as parameter 
a "vibrational" temperature T„ via: 



N 



dSv 
de 



(42) 



The calculation proceeds along the same line outlined in 
the case of the gaussian model. First, from the general- 
ization of Eq. ^2 i. e. from: 



1 
T 



dS/N dJ:{e)/N dSy/N 



de 



de 



de 



dT,{e)/N 1 

de 

(43) 

we get the temperature dependence of the energy of the 
visited minima: 



e(r) 



1 



1 



2 KksT ksTv 



(44) 



and inserting Eq. ^| into the definition of the gaussian 
model, Eq. 1131 we have the explicit expression of the 
configurational entropy as a function of the temperature: 



E(r) = ksN 



kfiT knT,, 



(45) 



We can now eliminate e by introducing the Kauzmann 
temperature defined by I](T;f ) = 0: 



e — 2^/akBTK 



Tv — Tk 



(46) 



thus, substituting this expression in Eq. 1451 

Tk 



E(T) = keNa 



1 



Tv — Tk 



T 



(47) 



Through the configurational entropy, we can apply 
Eq. ^|to find an expression for the fragility: 



TO, 



2TgTK{Tg/Tv) 



Tl)-2TK{Tg^TK){TjTv) 



(48) 



In this expression, besides Tg -the parameter that em- 
bodies our choice of the value of viscosity that define the 
glass transition temperature- there are the two system- 
dependent parameter Tk (a way to express e) and T^. Fi- 
nally the temperature dependence of the viscosity turns 
out to be controlled by the law: 



r/(T) = rjoo exp I 



DT 



K 



T 



T-TkT + Tk-2T{Tk/Tv) 

2 



Tk 



T, 



(49) 



with, as before, D = £ / {aNksTK)- 

Similarly to Fig. 13 in Fig. |^ we report the temperature 
dependence of the configurational entropy of the gaussian 
model with energy dependent vibrational entropy as a 
function of Tk/T for different values of T^/Tk (reported 
in the figure) and compared with the similar quantity for 
the gaussian and the hyperbolic models. 
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FIG. 6: Similar to Fig. \^ the temperature dependence of 
the configurational entropy is reported for the gaussian model 
with energy dependent vibrational entropy for different values 
of the parameter T^/Tk (dot-dashed line). For comparison, 
also the two corresponding function for the gaussian (full line) 
and hyperbolic (dashed line) models are reported. 

Analogously, Fig.Hshows the temperature dependence 
of the viscosity as predicted by the gaussian model with 
energy dependent vibrational entropy for different values 
of Tv/Tk- As can be noticed, it seems that the values 
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FIG. 7: Similar to Fig. |^ the temperature dependence of 
the viscosity is reported for the gaussian model with energy 
dependent vibrational entropy for different values of the pa- 
rameter Tv/Tk (dot-dashed line). For comparison, also the 
two corresponding function for the gaussian (full line) and 
hyperbolic (dashed line) models are reported. 



Speedy's expression of fragility 

In 1999 Speedy -working in the framework of the 
gaussian model and assuming the vahdity of the Adam- 
Gibbs relation- choose to express ("/" in his lan- 
guage) in term of a and '^{Tg) ("A^gSiTg)" in Ref. ,2^). 
With these variables, Eg. 1201 becomes: 



2a 



1. 



(50) 



Speedy used this relation to state that "..this quantifies 
the Angell observation that fragile liquids sample more 
basins in configuration space than strong liquids" . Actu- 
ally, Eq. I5UI docs not help much in establishing whether 
or not the reported Angell observation is correct. Indeed 
the proportionality between and a holds only if one 
neglects the possibility that S(Tg), a system-dependent 
quantity, depend on a. In principle its implicit depen- 
dence on a can also reverse the fragility-number of states 
relation. 



of allows one to interpolate between the behavior of 
the gaussian model (obviously reached for T^, — + cxd or 
— oo) and that of the hyperbolic model (that is approxi- 
mately obtained for T^/Tk ~ — 1 — 1.5). It is worth to 
remember that, in most numerical simulations of model 
liquids, is found to be negative for constant density 
(thus constant PEL) simulations, while T^, > for con- 
stant pressure simulations M- In the case of a model 
for water , the sig n of Ty has been found to be density de- 
pendent [43, 153I ■ On the experimental side, at constant 
pressure, the sign of turns out to be both positive [s^ 
and negative |53|, depending on the specific system. 

Finally in Table II we report the relevant expression 
relative to the gaussian model with energy dependent 
vibrational entropy {Ty ^ 00) compared with those of 
the gaussian model [T^ — 00). 



Sastry's expression of fragility 

More recently another expression for the fragilit y in 
term of the PEL features was derived by Sastry l26l|. 
Also in this case the gaussian model and the Adam-Gibbs 
equation are at the basis of the theory. However, Sastry 
does not use Eq. ^| to obtain the fragility. He assumed 
i) the validity of the VTF law, so to relate (compare 
Eq. El and ^ the configurational entropy to the coeffi- 
cient D, which, as discussed before, is an index of ki- 
netic fragility (actually, Sastry reports his expression for 
the fragility K = 1/-D), and ii) the coincidence of To 
with Tk. The Sastry expression takes also into account 
the possible energy-depth dependence of the basin vibra- 
tional free energy. In order to compare the expression 
reported in Ref. Ei^ with Eq.s 1201 and (201 however, we 
can put the quantity 5S (in Sastry's notation) equal to 
zero. The Sastry expression becomes (with the change of 
notation from cr to e): 



FRAGILITY AND NUMBER OF STATES 

In the following sections we will discuss the possibility 
to predict the fragility of a system from the knowledge 
of the parameters characterizing the distribution of the 
minima of the PEL. First, we analyze the recent works 
that have attempted to relate the fragility to the "num- 
ber of states". Secondly we will see how -given a fixed 
configurational entropy model- one can obtain the whole 
range of fragilities, thus demonstrating that, in order to 
asses the fragility of a system, some additional informa- 
tion is needed. 



K 



2£ 



1 



T„ 



(51) 



Here "Tg" is the MD glass transition temperatures. In 
Eqj^we have explicitly included the Adam-Gibbs con- 
stant £ which was implicitly a.ssumed constant and land- 
scape independent in Ref. ^^(see also [s^]). 

After the conversion from K to rrig , using Eq. [SI we 



have: 



TO, 



171n(10) 



To 



+ 1 



(52) 



Similarly to Eq. j^Ol also this equation cannot be used 
to predict the a dependence of the fragility. Indeed, a 
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ksTK 
e 

e{T) - eo 
•E{T)/NkB 



lnir,{T)/vo.) 



de 



= 



2^ 



2akBTK 



2k _ 



(Tg+T^ )-2Tg Tk (Tg / T„ ) 



T-Tfc T+Tj(-2T(Tj 



TABLE II: Summary of the main relations relating the relevant quantities (left column) for the gaussian configurational entropy 
models: the simple gaussian model {Tv — oo) and the gaussian model with energy dependent vibrational entropy. 



appears here explicitly but also implicitly, via the system- 
dependent quantities e and Tx (see Table I). Finally, we 
want to stress that the approach followed in the deriva- 
tion of the previous expression of the fragility is intrinsi- 
cally inconsistent. Indeed, as previously pointed out, the 
gaussian landscape (i), the VTF law (ii) and the Adam- 
Gibbs relation (iii) are not mutually consistent and, as 
also noticed by Sastry [i^ , the hypothesis i)-iii) can only 
be consistent if one uses a low temperature expansion of 

E(r). 



Can the fragility be derived entirely 
from the configurational entropy? 



We aim now to prove with an example that, in gen- 
eral, the configurational entropy alone is not sufficient 
to determine the fragility of a system. We will use the 
gaussian model for the configurational entropy and, with 
the help of Eq. Efl we will set-up an "AngeU plot". We 
could have selected any other landscape model, reach- 
ing the same conclusion. Let's suppose to have an hy- 
pothetical system, fully defined by a gaussian landscape 
with a given value of the relevant parameters a, e and 
Co- The temperature dependence of the viscosity in 
this model is reported in Eq. [2] To set up an AngeU 
plot, we need to define the "glass transition temperature" 
Tg. As done experimentally , once the T-dependence of 
the viscosity is known, Tg is defined from the condition 
log(?7(Tg)/77oo) = 17. Using Eq. EH the solution of this 



equation for (positive) Tg is: 



T 




171n(10) J 



(53) 



with D = S/iaNksTK) = 2£/{y/^N£). For sake of 
compactness, let us define the function j{x): 



-fix) 



so that 




2 171n(10) 



Tg = Tki{D). 



(54) 



(55) 



Obviously, the expression of Tg, besides the trivial tem- 
perature scale Tk, depends on the parameter D[= 
2£ / {^/aNe)) that, in turn, embodies the information on 
the " number of states" but also from quantities distinct 
from the statistic of the minima (specifically from the 
parameter £). We want now to plot the re-scaled loga- 
rithmic viscosity 2/(T) = \\og{r^{T)/'n^)]/[\og{r]{Tg)/'noo)] 
as a function of Tg/T . The quantity y{T), by definition 
ofTg, turns out to be equal to y{T) = [log(77(r)/77oo)]/17, 
or, by using the expression for ri{T) reported in Eq. 1211 
to: 

D Tk T 



y{T) 



(56) 



n \n{lQ) T -TkT + Tk 
We can now eliminate Tk from this equation in favor of 
Tg using Eq.ElKTK = Tgh{D)) to get: 



V{T) 



(57) 
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In Fig. I^lwe have reported the quantity y{T) of Eg. [571 
vs. Tg/T, i. e. we have made an AngeU plot, for different 
values of the parameter £ at fixed a and e. The fragili- 
ties TOs are the slopes of these curves at the upper right 
corner of the plot. What is remarkable here is that, by 
varying the quantity £ entering in the numerator of the 
exponent in the Adam-Gibbs relation (Eq.[9]) at fixed 
configurational entropy, we can span the whole range of 
fragilities. In other words, for a given (gaussian in the 
present example) landscape, with well defined statistical 
properties (fixed a and e), we can have a strong system 
(large £) as well as a fragile one (small £). Therefore, we 
conclude this section with the statement that in principle 
-whenever the Adam Gibbs relation represents a good ap- 
proximation of the relation between transport properties 
and configurational entropy- the knowledge of the config- 
urational entropy alone would be not sufficient to define 
the fragility of a system\56\. This statement, and the role 
of the effective barrier height in determining the fragility 
of a glass has been already discussed in literature (see 
e. g. Ref. ^57]). The previous conclusion does not imply 
that the fragility cannot be derived from the landscape 
properties: indeed, it is possible, and actually most likely, 
that the quantity £ could be derived from other features 
of the PEL than the minima distribution, as, for example 
the minimum-to-minimum barrier heights. Future stud- 
ies must focus on the relation between £ and the PEL 
properties and on the physical range of values of £. 




T /T 

a 
c 

FIG . 8 : Reconstructed Angell plot 

([log(r,(T)/r?oc)]/[log(r,(rg)/r,oo)l vs. TjT) for the case 
of the gaussian model. The different curves correspond to 
different £ values: from top to bottom £'=80, 40, 20, 10 and 
5 in units of ^/aNe/2. 



Strong-to-fragile transition 

In the previous section we have shown that, on a gen- 
eral ground, a simple gaussian landscape with fixed sta- 



tistical properties could be shared by the whole class of 
known systems; they would simply differ in the value of 
£ that, in turn, induces a different value of Tg/Tx, thus a 
different fragility. In this scheme a fragile system -having 
Tg close to Tk (as deduced from Ea. l^ - visits that part 
of the landscape where 5](e) is strongly e dependent, thus 
(see Eq. [Tn|) pushing up. On the contrary, a strong 
system has Tg far away from Tk, and the system is con- 
fined to visit the region where S(e) is almost flat. In 
other words, if all the system shared the same landscape, 
due to the difference in the parameter £, a strong system 
(large £) would visit the "top-of- the-landscape" , while 
a fragile system (small £) would be allowed to go down 
in energy. If this scenario were correct, we would ex- 
pect that real systems verify Eq. [201 (or [201 or I39|) . In 
Fig. [01 we report, for those systems where all the three 
quantities Tg, Tk and are known (see Table III), 
the fragility as a function of Tg/TK (symbols). Also 
shown in the same figure are the predictions of Eg. s 1201 
[201 [201 (lines). Few points must be underlined: i) there is 
a rather good general agreement, but the single systems 
does not strictly verifies none of the three predictions. 
This can be due to the existence of landscapes different 
from the three simple cases discussed at the beginning 
of this chapter, or, most likely, to the presence of a fi- 
nite value for Ty. Indeed, recent molecular simulations 
of model liquids 0| clearly show such a phenomenology, 
indicating a non negligible energy dependence of the vi- 
brational entropy, ii) The differences among Eq. 1201 and 
[201 are so small, that the experimental data do not al- 
low to discriminate among these two different landscape 
models, while the (pure) logarithmic model seems to be 
definitively unacceptable. Most likely, a gaussian model 
with a small logarithmic correction would be still accept- 
able, iii) Among the systems represented in Fig. [01 the 
lowest fragility is ~35, i. e. the strong systems are ab- 
sent (for these systems a reliable estimation of Tk does 
not exist), and this does not allow to firmly establish 
the general validity of one of the three model, and, more 
generally, of the idea presented before that strong sys- 
tems and fragile systems are characterized by a common 
configurational entropy and a different elevation in the 
PEL. 

Of course, we are not stating that the depicted be- 
haviour is the actual one. Different systems have different 
value of a and may even not be described by a common 
landscape model. A typical example, that it is worth to 
discuss here, is the case of silica. As shown by Horbach 
and Kob [^, vitreous silica (as described by the BKS 
potential model [s^) show strong T-dependence of the 
fragility. More specifically, v-Si02, which is a well known 
strong system close to the glass transition temperature, 
turns toward a more fragile behaviour on increasing T. 
This phenomenon, called " strong-to- fragile" transition, 
first proposed for the case of water |60j, has been observed 
in simulations of water |6l| and (simulated) berillium flu- 
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T, Tk Ref. 



2- mety Ip ent ane 


80.5 


58 
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[24], [24] 


butyronitrile 


100 


81.2 


47 


[241, [24] 


ethanol 


92.5 


71 


55 


[24],[55j 


n-propanol 


102.5 


73 


36.5 


[24], [24] 


toluene 
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96 
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[24],[55j 


1-2 propan diol 


172 


127 


52 
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glycerol 
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53 


[24], [24] 


triphenil phospate 
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160 


[24], [24] 


orthoterphenyl 
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81 


[2^, [2^ 


m-toluidine 


187 


154 


79 


[24], [24] 


propylene carbonate 


156 


127 
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[24], [24] 


sorbitol 


266 
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93 


[24],[24] 


selenium 


oU ( 


z4U 


1 


[2^,|i4{ 


ZnCl2 


380 


250 


30 


[24],[55j 


AS2S3 


455 


265 


36 


[24], [55] 


CaAlaSizOg 


1118 


815 


53 


mm 


Propilen Glycol 


167 


127 


52 


[27],[27j 


3-Methyl pentane 


77 


58.4 


36 


[27], [27] 


3-Bromopentane 


108 


82.5 


53 


[27], [22 


2-methyltetrahydrofuran 


91 


69.3 


65 


[27], [27] 



TABLE III: Summary of the quantity relevant to test Eg. 1201 iTr.. Tk and rUg) for those systems where the three quantity are 
all known. The last column reports the references where the data has been found. In those case where more than one value of 
the parameters are known, we have reported here the average value. The first of the two references refers to the couple (Tg, 
Tk) and the second to m^. 



oride [62| ■ It is obvious that the fragile-to-strong transi- 
tion cannot be framed within the possibility described in 
the first paragraph of this section. 

In a recent simulations work, Saika-Voivod, Sciortino 
and Poole [H^l have shown that the configurational en- 
tropy for liquid sihca -as derived from a MD simula- 
tion based on the BKS interaction potential model- is 
far from being " gaussian" . More specificaUy, they found 
that S(r) -at low T- shows a tendency towards a posi- 
tive curvature and does not seem to extrapolate to zero 
entropy at a finite temperature. This behavior is shared 
by the logarithmic model (or by a combination of the log- 
arithmic and gaussian models, as proposed in Ref. |49|) 
for S(e). This model predicts an infinite slope of E(e) 
at ck, and this could be in agreement with the simula- 
tion results of Ref. [g^ as the low statistics in the tail 
of Vl]si{E) -as measured by MD- does not allow to safely 
determine dT,(e)/de evaluated at ex- 

The logarithmic model, however, similarly to the other 
models presented before, is not capable to catch the 
physics of the strong-to-fragile transition. Indeed, the 
fragility expressions for all the examined models fEa.s l20l 
1301 139|) show a monotonic T dependence, whit a tendency 



toward a decrease of the fragility on increasing tempera- 
ture (see Fig. I10|l . A behavior opposite to what observed 
in simulated vitreous silica. It is therefore clear that 
an infinite slope of S(e) at bk alone is not sufficient to 
guarantee the existence of a strong-to-fragile transition. 
What is actually sufficient (necessary ?) for a strong to 
fragile transition -i. e. to have a maximum in the vs. 
Tg/Tji function- is that the configurational entropy -as a 
function of T- had a non-zero limit for T 0. This can 
be understood looking at Eq. ^| It is clear that a fragile 
system is characterized by a large value of I]'(T) (fragile 
systems explore the "steep" part of the PEL), while a 
strong system will have a small value of S'(T), but also a 
non zero S(r). This certainly happens at the "top of the 
landscape", but could also happen at low T if S(0) 7^ 
(in the logarithmic model, at low T, Y^'{T) 0, but 
the same does S(T) and the resulting fragility increases 
continuously). Thus, a strong-to- fragile transition could 
take place only if the landscape of the systems allows for 
a finite number of states at zero temperature, i. e. for 
a (exponentially large with N) degenerate fundamental 
state. The existence of such a degeneracy for system with 
short range interaction (non mean field systems) poses 
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FIG. 9: Experimental values of the kinetic fragility plot- 
ted as a function of the ratio T^/Tk for those systems where 
the three quantities (m^, Tg, and Tk) are available. The in- 
put data are reported in Table III. For those systems where 
more than one determination of the parameters is known, we 
have reported in the plot the average value together with an 
"error" bar that indicates the whole dispersion. 

several problems (see the discussion in Ref. 0), and is 
certainly calling for further investigation. 
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FIG. 10; Temperature dependence of the fragility for the three 
examined models: gaussian (full line), hyperbolic (dashed 
line) and logarithmic (dot-dashed line). The quantity is 
reported as a function of Tg/TK- 



DISCUSSION AND CONCLUSION 

In conclusion, in this paper we have first summarized 
the main definition of fragility, then we have recalled and 
studied different models for the configurational entropies 
present in the literature. Using the Adam-Gibbs relation 
to link the dynamics of a glass forming system to its con- 
figurational entropy, we have reported the explicit expres- 



sions for different quantities, among which the fragility. 
From the reported relation, it is clear that in general 
the fragility cannot he derived by the knowledge of the 
configurational entropy. More specifically, given a fixed 
" landscape" , different system fragility can be mimicked 
by varying the parameter £ entering in the numerator of 
the exponent of the Adam-Gibbs equation. On a general 
ground, the fragility of a system depends on the ratio 
S/aNksTK. 

The fact that the whole range of fragility can be de- 
rived from a given PEL model (e. g. the gaussian model) 
with the same statistical properties seems an interesting 
possibility. If this was the case, the strong glass forming 
materials would be characterized by a large value of £ and 
would explore the "top-of-the-landscape" , while the most 
fragile ones would have small £ and would visit the states 
around the inflection point of S(T). Obviously other pos- 
sibilities exist, as for example that all the systems were 
characterized by the same £, and in this case strong glass 
would have a small number of states (small a), at vari- 
ance to the fragile systems with more states (large a). 
A further scenario can be hypothesized, that would also 
explain the existence of a strong-to-fragile transition: in 
this case the strong systems would explore the bottom of 
a landscape characterized by a non- vanishing zero point 
entropy. This is an interesting possibility that deserve 
deeper investigation. 

Overall, the present discussion, which heavily build on 
the validity of the Adam and Gibbs relation, indicates 
that in principle at least two possible classes of strong 
glass forming materials can actually exist. On one side 
we have those systems that -close to Tg- visit state at 
the top of the landscape and have a " regular" (gaussian- 
like) configurational entropy (let's call these systems as 
class A strong glass forming materials). On the other 
side we find the -let's say- class B strong liquids, that 
visit minima deep in the PEL, but with an exponentially 
large degeneracy of the fundamental state. The answer 
to the question whether class A and/or B strong systems 
actually exist requires further investigations. 

As a final comment, we would like to recall that 
fragility is often measured at constant pressure, while all 
the configurational entropy based models -as those pre- 
sented here- are built on the assumption of a well defined 
S(e) function, i. e. they assume constant density. The 
relationship between constant density and constant pres- 
sure fragilities is one of the topic under discussion at the 
present time. As an exam ple , in the case of soft sphere 
systems it has been shown that T,c{T) along isochoric 
and isobaric paths are very close one to each other. Sim- 
ilarly, in a very recent work 65], Tarjus and coworkers 
show that -in alcohols- the change in density only slightly 
affects the fragility, thus indicating that under the exper- 
imentally accessible density changes the landscape suffers 
only minor modifications. For other systems, on the con- 
trary, it has been observed a large deviation (« 40%) 
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between constant density and constant pressure fragili- 
ties 01 . This ongoing discussions, however, does not 
affect the conclusions of the present work since all the 
formalism could have been based on the enthalpy land- 
scape, instead of energy landscape, without any changes 
in the results. 
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